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<N 

The Falicov-Kimball model is a lattice model of itinerant spinless fermions ("electrons") inter- 
acting by an on-site potential with classical particles ("ions"). We continue the investigations of the 
I crystalline ground states that appear for various filling of electrons and ions, for large coupling. We 

investigate the model for square as well as triangular lattices. New ground states are found and the 
effects of a magnetic fiux on the structure of the phase diagram is studied. The fiux phase problem 



where one has to find the optimal fiux configurations and the nuclei configurations is also solved 



I in some cases. Finaly we consider a model where the fermions are replaced by hard-core bosons. 

^ This model also has crystalline ground states. Therefore their existence does not require the Pauli 

O principle, but only the on-site hard-core constraint for the itinerant particles. 

o 
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1 Introduction 



The occurence of long-range order in itinerant electron models, like the Hubbard 
model, is usually very difficult to prove. At half-filling on a cubic two or three dimen- 
sional lattice the Hubbard model is supposed to display Neel order at low temperature 
[d > 3) or in the ground state [d = 2). A simpler model where a proof of this fact has 
been achieved is the so called Falicov-Kimball [12]. There, only one species of elec- 
trons (say the spin up electrons) is allowed to hop, while the other type of electrons 
(spin down electrons) are treated as a classical lattice gas. A nice interpretation of the 
model is that of a model for matter, where the static particles correspond to heavy 
ions which interact with the electrons via an on-site potential. The system displays 
crystalline long range order for the ions which comes from the interplay between the 
on-site interaction of electrons and ions and the electronic kinetic energy. The Fermi 
statistics of the electrons plays an important role, because if the electrons are replaced 
by regular bosons then the ions prefer to clump together [12]. An open question which 
we solve in this work is what happens if fermions are replaced by hard-core bosons. 

On a bipartite lattice and for a special value of the chemical potentials (/ie, f^i) of 
the electrons and ions, it was shown in [12] that at low temperatures the ions tend to 
occupy one of the two sublattices. On a square lattice this corresponds to a chessboard 
configuration. This value of (/ie,/Wi) is such that the system possesses a particle-hole 
symmetry. In what follows we fix this special value at (/ie,/Wi) = (0,0) and refer to 
it as the symmetry point (see paragraph 2.1). This result was then extended to a 
narrow strip of chemical potentials around the symmetry point, provided the coupling 
between electrons and nuclei is strong [16]. 

Concerning the ground states of the model more results are available. The study 
of the zero temperature phase diagram in the chemical potential plane (/ie,/Wi) was 
undertook in [7, 8] and for one dimension in [9]. For large coupling, which is the case 
that will interest us in this paper, formal perturbation expansion shows that away from 
the symmetry point new periodic ground states appear, with rational densities of ions 
different from one half. The existence of some of these ground states was confirmed 
rigorously in [11], where a canonical setting was used. For the fixed densities 1/3, 1/4, 
1/5, the ground state configurations are determined rigorously using perturbation 
theory for large coupling. 

In the present work we continue the investigation of the phase diagram in the 
chemical potential plane (/ie,/Wi) in a more systematic way. Our results generalize 
and extend what was previously known in three directions : i) we consider square as 
well as triangular lattices, ii) we allow for Fermi statistics for the electrons as well as 
hard-core boson statistics, iii) we investigate the effects of a magnetic field which can 
be homogeneous or inhomogeneous. 

Let us briefiy describe the new features which appear. 

i) Triangular lattice. To our knowledge nothing rigorous is known for the Falicov- 
Kimball model on a triangular lattice. This lattice is not bipartite and in general the 
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particle-hole symmetry is not present. However it is recovered for a special magnetic 
field, namely a flux equal to ±^ through each triangle (see the end of paragraph 2.1 
for more details). This fact has been exploited recently in [24] to extend the uniform 
density theorem [20] to non bipartite lattices with special magnetic fluxes. To leading 
order of the perturbation theory one sees that the effective interaction between the 
nuclei is of the Ising type (for any lattice). Thus when the lattice is triangular there 
appears frustration effects at zero chemical potential /ie = f^i = 0. An interesting 
question is therefore what happens when the higher order terms of the perturbation 
theory are included : does this frustration disappears or does it persist ? Our study 
indicates that frustration disappear when higher terms are included. In fact for zero 
magnetic field it is rigorously established that there is no frustration at /ie = fii = 
for large enough coupling constant. 

ii) Quantum statistics. The "electrons" will either be regular fermionic particles 
or hard-core bosons. For usual bosons, as stated before, there can be no ordering 
for the ions which prefer to clump together. On the contrary, for hard-core bosons 
there appears crystalline order for a finite range of densities. The ground states may 
be different from those of fermions. However for the square lattice around the sym- 
metry point the situation is the same as with the Fermi statistics : the ground state 
configuration of the ions are the chessboard states. Here the methods that we use are 
restricted to large coupling, but this last result can be proved to hold for all couplings, 
by using a refiection positivity technique [22]. 

In the case of Fermi statistics, given a configuration of ions, we can always reduce 
the problem to a one particle hamiltonian. One can then compute the ground state 
energy simply by filling the levels (they depend on the configuration of nuclei) of the 
one particle hamiltonian according to the Pauli principle. For hard-core bosons the 
situation is completely different a priori. The hard-core introduces a non trivial two 
body interaction between the bosons. Therefore we have a genuinely many particle 
hamiltonian. To perform the perturbation expansion we use a functional representa- 
tion recently developped by S. Miracle-Sole and A. Messager [25] for the model in the 
fermionic case. In fact for hard-core bosons the method is perfectly well adapted and is 
even simpler. It turns out that the effective potential between the ions is qualitatively 
similar to that of the fermionic case, and also leads to crystalline ground states. 

The chemical potentials plane is decomposed into domains corresponding to a 
fixed density of electrons and ions. Therefore the compressibility vanishes at zero 
temperature and these states are insulating. It is interesting to note that in the 
bosonic case this behaviour is proven for a truly many particle interacting system. 

iii) Magnetic fields. The hopping amplitude of the itinerant fermions or bosons can 
be complex. The sum of the phases around a circuit of the lattice can be interpreted 
as the fiux of an external magnetic field. Our results allow to see how the ground 
states are affected when one varies the magnetic field. It is found that in general the 
crystalline configurations of the nuclei can vary with the magnetic field. In fact this is 
not surprising since in this model the crystalline order is produced by a combination 
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of the potential and the kinetic energies. 

Another problem that we address is the so-called flux phase problem. This prob- 
lem arises in the context of theories for high Tc superconductors [f, f5], but is also 
important from the point of view of understanding electronic diamagnetism of itiner- 
ant electrons on two-dimensional lattices [fO]. In the present context one fixes the fiux 
configuration through each face (plaquette or triangle) of the lattice. The problem is 
to minimize the ground state energy simultaneously with respect to the ion and fiux 
configurations. In the bosonic case it is easy to see that the optimal fiux is always zero 
(diamagnetic inequality). For fermions however the situation is more subtle. We find 
that for large coupling, the optimal fiux configuration consists of zero and tt through 
the faces of the lattice. On the square lattice, for density | the optimal fiux configu- 
ration is TT through each plaquette; for density | the optimal fiux configuration is not 
uniform, but of period three, the average fiux being ^; for density | one finds that the 
fiux is for sufficiently large coupling. We refer the reader to section 4 for more details 
and the results on a triangular lattice. It has been conjectured in the literature that 
for free fermions on a square lattice the optimal fiux should be equal to 27rn, where 
< n < 1 is the filling factor [10]. Our results cannot however be compared with this 
conjecture since we work at large coupling, and we are therefore far from the limit of 
free electrons. On the other hand they show that the optimal fiux configuration is not 
necessarily uniform as it is often assumed in the literature. Recently the conjecture 
has been proven for the half-filled band of the Hubbard model [18] (see also [23]). 

The paper is organized as follows. In section 2 we derive the expansions 
for the fermionic and bosonic cases for square and triangular lattices. The phase 
diagrams for a uniform magnetic fiux are then studied in section 3. Section 4 discusses 
the optimal fiux configurations of the model. The appendices contain more technical 
material. 

2 — expansion for the Falicov— Kimball 

Models 

2.1 Definition of the models 

The spinless Falicov-Kimball models have been introduced to discuss lattice sys- 
tems consisting of identical quantum particles interacting with classical particles (or 
"ions") by means of on-site interaction which mimics the Coulomb repulsion or at- 
traction. The quantum particles can be subjected to an external magnetic field which 
appears in their kinetic energy. We shall investigate several models in which the quan- 
tum particles can be either fermions or hard-core bosons; the lattice on which the 
systems are defined is either a square or a triangular lattice. 

We shall consider only finite lattices A and for technical reasons we introduce 
periodic boundary conditions. The quantum particles are represented by creation and 
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annihilation operators, which satisfy the usual "commutation" relations : 

fermions : Wt ^ ^ij'} = Wx, dy} = Vx, j/ G A 

{a^,a+} = S^y (2.1) 
hard-core bosons : [a*, a*] = Vx ^ y ^ A 

{atf = {a,f = {«.,«+} = ! (2.2) 

The "ions" are described by a random variable G { — 1,+!} such that s^; = +1 if 
the site x G A is occupied and s^; = — 1 if the site x is empty. The finite volume 
hamiltonian of the Falicov-Kimball models is written as 

Ha{s) = - J2 ^xy^tay + UJ2 ^Aata^ - ^) (2.3) 

where s = {sx}xeA denotes the configuration of ions, t^y = \txy\e^^''^ = t*^, i.e. O^y = 
—Oyx^ and 6xy G [0, 27r[ is related to the integral of the vector potential O^y = /J ^ ■ dx. 
For the cases under investigation \txy \ = t for nearest neighbours and zero otherwise. 

Given a circuit io of length |c<;| = n, i.e. an ordered sequence of sites (xi, . . . , a^i), 
with txjXj+i 7^ for j = 1, . . . , n, defined up to a cyclic permutation, we introduce the 
fiux across the circuit as 

n 

<^c. = E^-.-.+i (mod27r). (2.4) 
1=1 

Moreover, we also need to introduce the number N[ijj) of periodic sequences in cu, i.e. 

Niu) = ^ (2.5) 
n' 

where n' is the length of the shortest circuit io' such that io is obtained by repeating 
A^(c<;)-times the circuit io'. 

If the coupling U is positive the interaction between quantum particles and ions 
is repulsive, if U is negative it is attractive. 

It is important to stress that contrary to the case of fermions where, except for 
the statistics, we have a system of non-interacting particles under the infiuence of the 
external potential {Usx} and the magnetic field, in the case of bosons we have a system 
of interacting particles. Therefore in the case of fermions one is led to consider the 
eigenvalues of the one particle hamiltonian, but in the case of bosons we are forced to 
consider the N particle hamiltonian. 

As usual, the effective interaction between the ions is defined by the trace of 
the Boltzmann factor over the Fock space of the quantum particles 

g-/3FA(/3,s;Ate,A'0 _ 'Pj,|g-/?[-ffA(s)-AteiVe-A'«JV«(s)] | (^2 6) 

where 

iVe = ^a+a. N,{s) = W{s^ + 1) (2.7) 
xeA xeA 
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are the electron and ion numbers, and /ie, f^i their respective chemical potentials. The 
zero-temperature limit defines the effective hamiltonian for the ions 

Heff{s; fie, fii) = lim Fa{/3; s; fi^, /«0- (2-8) 

p^oo 

In the case of fermions the particles do not interact and the hamiltonian (2.3) is 
the second-quantized version of the one-body hamiltonian (up to the constant term) 

hA{s) = T + US (2.9) 

{T^xy ^xy {^)xy ^x^xy (2.10) 

therefore the effective hamiltonian (2.8) for the ions is the sum of the lowest eigenvalues 
ej(s) of h^is) 

Heffis; /ie, Hi) = J2 i^ii^) - /^e) - -^TrS* - fi,N,{s) 

ej(s)</ie 

= E{s- fie) - fieNe - iJ.,N,{s) (2.11) 

with 

E{s;fie)= J2 e,(3)-|Tr^, iV^ = ^ 1. (2.12) 

ej(s)</ie ej(s)</ie 

Let $ = {(f>c} denote the fiux configuration across elementary cells C (i.e. the pla- 
quettes P or triangles A of A). When (pc = (p for all C we say that $ is homogeneous. 
The problems we want to investigate are the following. 

i) given (/ie, /Wi, $), with $ homogeneous, find the configurations of ions S = {s"} 

which minimize the effective hamiltonian : this is the problem of the ground 
state, with constant external magnetic field. 

ii) given (/ie,/Wi), find the configurations of ions S and the fiuxes {^c}, not neces- 

sarily uniform, which minimize the effective hamiltonian : this is the fiux phase 
problem. 

We shall use the expansion of the effective hamiltonian Hefj and we shall 
restrict our considerations to the case 

\U\>zt fie(E]-\U\ + zt,\U\-zt[, t = sup\txy\ (2.13) 

where z is the maximal coordination number [z = A for the square lattice, z = 6 
for the triangular lattice). This is the domain where the expansion is absolutely and 
uniformly convergent with respect to A. 

Moreover, if U > (repulsion), (2.13) implies the ^^half-filling condition^^ Ne + 
Ni[s) = |A| and if U < (attraction), (2.13) implies the neutrality condition^^ = 
NAs). 
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This is a direct consequence of the properties of the spectrum of the one-body 
hamiltonian h^is) (for fermions) and also follows from the closed loop expansion (sec- 
tion 2.3) valid for fermions and bosons. 

Finally we want to stress the role played by the statistics and the nature of the 
lattice on the symmetries of the phase diagrams. First of all if s = {sx}xeA is a ground 
state for U and (/ie,/Wi), then s = {—Sx}xeA is a ground state for —U and (/ie,— /Wi). 
This is true for any lattice, any statistics and any flux configuration. Therefore we 
can restrict our discussion to the case [/ > in the next sections. 

For bipartite lattices and fermions (the square lattice for example) the phase 
diagram has the following symmetry. If s is a ground state for some value (/ie,/Wi) 
and some fiux configuration {^p}, then s is a ground state for (— /Ue, and { — (pp}. 
This follows from the unitary particle-hole transformation a'^ ^ a^^ ^ applied 
to the hamiltonian (2.3). By a time reversal transformation (complex conjugation) we 
see that s is also a ground state for (— /ie, — /Wi) and {^p}. In particular this implies 
that figures 7 and 15 in sections 3 and 4 are symmetric with respect to the line // = 0. 

For the triangular lattice and fermions (non bipartite case) the above symmetry 
of the phase diagram is in general lost. However the following holds. If s is a ground 
state for some value (/ie,/Wi) and some fiux configuration {^a}, then s is a ground 
state for (— /ie,— /Wi) and {tt — ^a}- This follows again from the same particle-hole 
transformation. By time reversal the statement still holds with {vr — ^a} replaced 
by {tt + ^a}- In particular if ^a = if through all triangles then tt — (pA = if, 
and we have a symmetry of the hamiltonian as mentioned in the introduction. As 
a consequence on the figures 10, 11, 12 the horizontal line ^ = f is symmetric with 
respect to the vertical line // = 0. In figure 16 the case // > can be obtained from 
/i < by interchanging the values and tt for the fiuxes and replacing s by s. 

In the case of hard-core bosons the particle-hole transformation a'^ ^ a^^ ^ 
leads to the following property, irrespective of the type of lattice. If s is a ground 
state for (/ie,/Wi) and {^c} then s is a ground state for (— /ie, — /Wi) and { — ^c*}. By 
time reversal s is also a ground state for (— /Ue, "/"O and {^c}- Therefore the phase 
diagrams of figures 8 and 13 are symmetric around the line // = 0. 



2.2 f7 ^ expansion of -ffe// for fermions 

For fermions systems satisfying the condition (2.13), we can use the representation 

E{s- /ie) = ^ f dzT, [ - ^Tr^A (2.14) 

27rz Jv \z — ht\s) j 2 

where I) is a circle in the complex 2;-plane, centered at and enclosing all negative 

eigenvalues. By expanding the resolvent we obtain the following absolutely convergent 
expansion in powers of \J~^ 
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kimr\i^J ^ '--•••'--(mH-i)!(^-mH-i)!ivH ^'-''^ 

where 

r |A| - N,{s) iiU>zt 
' \ N,{s) iiU < -zt. 

In the above formula ^k{s) = {cu G ^{s)] \lo\ = A;}, where 0(s) is the set of circuits 
Lo = (xi, . . . , Xfc, Xi) which contain at least one empty site and one occupied site. For 
[/ > 0, m[ijj) is the number of sites in (xi, . . . , x^) with s^^ = — 1, i.e. the number of 
empty sites in io; for [/ < 0, m[ijj) is the number of occupied sites in (xi, . . . , Xk). 

In the particular case where \txy\ = t if (x^y) are nearest neighbours and zero 
otherwise (2.15) yields 

+ Y.T7JTi-I gcois) cos ^Js) (2.16) 

k>2 1^ I LoeQkis) 

where 

a (s) = (-lY^i^)— {_k^^2)l 1_ 

^"^^ ^ ' 21^-^ {m{Lo)-iy.{k-m{Lo)-iy.N{Lo) ^ ' 

To conclude this discussion let us note that for bipartite lattices the sum in Eq. 
(2.15) is only over even values of k. 

2.3 Closed loop expansion for F 

For hard-core bosons a similar expansion is valid, although its derivation 
is very different. We use a method first invented (for fermions) by A. Messager and 
S. Miracle-Sole [25]. The basic steps are illustrated here for the convenience of the 
reader, in the case [/ > 0. 

Thanks to the Trotter product formula 



M^oo 



M ' 



f2.18) 



with Hq the kinetic part of (2.3) and Hi the interaction part. We can express the 
trace over the Fock space of the quantum particles in the form 



iVe>0 

M 

• lim ^ £e-^^^E^:^^(^')n«(^.;^.-i) (2.19) 



where 
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X = {xi, . . . , xn^} = ordered subset of A 
v[X) = \{x ^ X; Sx = l}\ = number of doubly occupied sites 

(by one ion and one electron) 
Xm = Xq up to permutations 

r 1 ifY = x 

a.[Y; X) = I '^tyx if Y differs from X by one jump from x io y 

I otherwise 
e = +f for hard-core bosons 

e = sign of the permutation Xq — )• Xm for fermions. 

Introducing the (I + z/)-dimensional lattice, where the new dimension is related to 
the "time" variable = -^n, n = 0. . . . , Af — I, and all the functions defined on this 
lattice are periodic in x^ with period /3, one is led to consider the family of oriented 
closed loops on the cylinder [0, Af — I] X A (not to be confused with the circuits defined 
previously). The following expansion can be obtained for \txy \ = t 

F{/3-s-iJ.,,iJ.i) = (/ie-/iOiV. -(/ie + |)|A| + F'(/3;3;/ie) (2.20) 
m^;/^e) = -^^im logl^* n KicoA (2.21) 

A^(cu,) = e-M(f^-A'e)KTle-^(^+^«)Kil |I e''^("^)e(cu,) (2.22) 

The sum * in (2.21) indicates the sum over non intersecting set of closed loops 
CO = (cui, . . . ^iOr) on the (I + z/)-dimensional lattice with at most one jump at each 
time Moreover ~ s means that the closed loops (cui , . . . , cu^ ) are compatible 
with s, i.e. the orientation of the segments [(x°, x), (a;° + x)] is in the positive time 
direction if s^; = +1, in the negative time direction if s^; = —I; and denote 

the number of segments in the time direction, which are oriented in the positive and 
in the negative time direction : X]j[|^jtI ~ l^jil] — \M ~ — hi^Uj) is the number 
of horizontal bonds (or jumps) in uf^ ^(^) = Zlj dxjXjj^i is the fiux through the circuit 
), projection of on A. 

We recognize in (2.21) the partition function of a polymer system in (I + p) 
dimensions. The polymers are the loops cOj and they cannot intersect. 

Using the standard cluster expansion we obtain the expansion of the effective 
interaction F in terms of connected closed loops 

OO 1 

-/3F'(/3;3;/ie)= lim ^- ^ [co, , . . . , co.,) (2.23) 

r = l UJi,...,UJr 

where is the truncated function. It has been shown in [25] that the above cluster 
expansion is absolutely convergent, uniformly with respect to M and /3, if [/ > Azt 
and \fie\ < U — 4:zt. 
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Remarks 

1. In one dimension e(cu) = 1 for all u and therefore we recover the known result 
that hard-core bosons are identical to fermions in one dimension (with open 
boundary conditions). 

2. In the case of fermions the condition of convergence |/ie| < U — zt, implies that 
/ie is in the gap separating positive and negative eigenvalues of hf^{s). 

3. In the zero temperature limit, and |/ie| < U — Azt^ (2.21) implies 

Heff{s; fj^e, IJ.i) = -{jli - /ie)iVj(s) - {jl^ + ^)|A| 




(2.24) 



and only loops which do not go around the cylinder [0, Af — 1] X A have to be 
considered. Indeed for any loop going around the cylinder there is a factor ^e~^^ 
which gives zero in the limit /3 — )• oo. 

Therefore in the limit /3 — )• oo, for each closed loop we have \ = and 
Ne + Ni[s) = |A| which is the half-filling condition (as it is expected since /ie is 
in the gap). 

2.4 Flux phase problem for bosons 

In (2.23), the right hand side is increased if each term of the sum is replaced 
by its modulus. This amounts to set e = +1 and to replace t^y by \txy\^ i.e. to set 
Gxy = 0. Therefore exp(— /3F') is increased if we replace the fiux by (pc = 0. The 
effective interaction, and as a consequence the ground state energy, are minimized if 
all the fiuxes are zero, for all values of (/ie, /Wi, fi). This property of bosons is in fact 
very general (see for example [26]). For fermions the situation is more subtle and is 
the subject of section 4. 

2.5 Expansion up to order 3 for fermions systems 

In the case of fermions systems, using Eq. (2.15) one obtains immediately the 
first few terms of the effective hamiltonian of ions, for U positive. Here the fiux 
configuration is uniform. 



Square lattice 



Heffis;^e,^i) = -7:ii^^ " ^e) ^ - ^(/ie + ^^ + U)\A\ 

^ xeA 



+ E 

<x,y>CA 



I^-tI^(^ + 2cos, 



+ E 



(x,y)cA 
\x-y\=V2 



16[/3 



(4 - cos (f))s^Sy 
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\x-y\ = 2 



(2.25) 



where the last sum runs over the plaquettes P consisting ol lour sites, and sp is the 
product ol the lour spins around the plaquette. 

Triangular lattice 

1 3 1 



+ E 

<x,y>CA 



■Sx-Sy + 2^ _|^g^3 



(4 — COS 2(f>)sxS^ 



(x,y)cA 
\x-y\=V3 



+ E 



(x,y)cA ACA ° 

k-!/|=2 



(2.26) 



with A a triangle and P two triangles sharing an edge; sa, -sp are the product ol their 
respective spins. 



2.6 Expansion up to order 3 for bosons systems 

For bosons systems we have to use the closed loop expansion (2.24). Let us first 
note that using the remark 3 ol Sec. 2.3, the activity (2.22) ol the loop io is 



(2.27) 



We then observe Irom Eq. (2.23) that the contribution in will come Irom 
loops with (k + 1) horizontal jumps (i.e. the coetficient ol t^~^^). 

Let us denote by (xi, . . . , Xk) the circuit defined by the projection (along the time 
direction) ol the loops (cui, . . . ,0;^) on A. In this manner the sum over the loops with 
(k + 1) jumps can be expressed as a sum over circuits. Now il lor a given circuit 
(xi, . . . , Xk) all the possible lamilies ol closed loops (cui, . . . ,0;^) associated with this 
circuit have the same lactor e(cu) = 11^=1 ^(^j)? then the contribution ol those loops 
lor bosons systems is immediately obtained Irom those already obtained lor lermions 
systems. 

In particular il the circuit (xi, . . . , Xk) is "linear" (does not enclose a plaquette), 
then one can check that e(cu) = 1. Furthermore il the circuit is sell-avoiding [xj ^ Xk 
lor all j ^ k), then there is exactly one loop and in this case e(cu) = e™~^ (e = +1 
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for bosons, e = — 1 for fermions). Therefore, in this case e(cu) = 1 for bosons and 
£{io) = — for fermions. 

These situations will hold for the expansion up to order 3 and therefore we obtain 
immediately the following expressions for bosons systems. One should note that the 
above simple properties which enable us to obtain directly the expansion for bosons 
from the corresponding expression for fermions do not work any more already at order 
4 and one has then to use Eq. (2.24). 

Square lattice 



Heff{s;fle,Hi) = -^(/^. - He) ^ 5^ - ^(/ie + /i, + U)\A\ 



+ E 



iU 16[/3 



(7 - 2 cos 



+ E 



{x,y)cA 
\x-y\=V2 



16[/3 



(4 + cos (f))s^Sy 



+ E 

(x,y)cA 
\x-y\ = 2 



f2.28) 



Triangular lattice 



1 1 3 ^3 

Heff{s;iJ^e,IJ'i) = --(/^. - /^e) E '^^ " g^^^' + P'^ + U + ^Tp)\M 

xEA 



+ E 

<x,y>CA 

+ E 



TT7 + 



cos > 



(7 — 5 cos 2^ 



^ (4 + cos 2(/))s^s,, + — rs^s,, 

(x,y)CA (x,y)CA 
\x-y\=V3 \x-y\ = 2 



(2.29) 



3 Phase diagram 
3.1 Method 

In this section we investigate the zero temperature phase diagram, i.e. the ground 
states of the effective hamiltonian, assuming a constant flux across each plaquette or 
triangle. In the next section we shall then solve the flux problem, namely given (/ie, f^i) 
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we shall find the set of fiuxes, and the configurations, corresponding to the minimum 
energy. 

We restrict our considerations to the case U > Azt and |/Ue| < U — Azt for which the 
previous formulae are valid. In this case the system satisfies the half-filling condition 
Ni-\-Ne = |A|, the expansion converges, and the effective hamiltonian is a function 
of [fii — fie) only. We then take /ie = and write fii = fi (which is to say fi = fii — fie)- 

The idea is to decompose Hefj at order k into two parts, the truncated hamil- 
tonian H^'^'^ (which depends on //), and the rest R^'^'^ (which is independent of //), 
i.e. 

Heff = H^'^ + R^'l 



H''' = -f E(^^ + 1) + E 777-- R'"' = E 777-. (3-1) 

^ X g>i q = k+l ^ 

In (3.1) we have dropped the constant — |[/|A|. Following [11] the truncated part can 
be written as H^^"^ = J2b Hb = ZIb H'^ where Hb is given by the expansion, 
5 C A with diam B < \{k + 1), and H'g = Hb + Kb- The "zero-potentials" Kb 
satisfy 

Y,Kb{s) = for alls (3.2) 

B 

and have the symmetries of the lattice. 

We shall then show that there exist fii < fi2 < - - - < fin (n- and fij depending on 
the order k and the fiux configuration $) such that for the following results hold : 

for fi e]fij,fij+i[, j = 0,. . . ,n, fio = -oo, fi^+i = +oo, 

a) it is possible to find Kb (depending on fi and $) such that H'^ is an m-potential, 
i.e. there exists a set of periodic configurations Sj = {s"}, where the s" are 
related by the symmetries of the lattice, satisfying the following conditions 

f H'sis) = H'sisJ) = min, H'sirs) if s\b = -sJIb ^3 3^ 



H'b{s) > H'sisJ) H-sIb^-sJ 



\B 



Therefore if fi G]/Uj, fij+i[^ the set of configurations Sj is the set of ground states 
for the truncated hamiltonian H^'^K On the other hand, for // = fi.j the number 
of ground states for H^'^'^ is infinite. 

b) it is possible to find e^, 6j positive, -|- 6j < fij+i — fij^ such that for fi G 
]fij -\- ej,/ij_|_i — (5j[, Sj is the set of ground states for the effective hamiltonian 
Heff- To establish this result, we follow the ideas introduced in [11] and write 
the rest as 

= (3.4) 
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where 

and R^J''^ is the average value ol R^J''^ over a cell C'p[x)^ centered at x, associated 
with the translation symmetry group of the s" (the precise definition of Cp{x) is 
given later), i.e. 

It is then easy to see that 

R^Hsf) = B^Hsr), for all 3f G <S,. (3.7) 
c) the /Uj, J = 1, 2, . . . , n, are functions of U and k of the form 

k ^q + 1 
g=0 ^ 

and Sj, Sj are of the order t^~^'^ / U^~^^ . 
To illustrate the method we consider the order : 

Hen = -f E(-^ + 1) + E = + (3-9) 

The ground states of H^'^^ are evidently 

S- = where s_ = {s^ = —1 Vx} for // < 

iS+ = where s+ = {s^; = +1 Vx} for // > 

For the domain // < 0, let // = —(5 and T[s) = {x : = +1}, 

i7e//(3) - H^ff{s_) = S\T{s)\ + E 777 hi^) - 

g>l 



>|r(. 



g>l ^ 



(3.10) 



Therefore for // < — 2;|^(1 — z-^) ^, s_ is the only ground state for Hpfj. Similarly for 
fi > zjj{l — zjj)~^, s+ is the only ground state for Hefj. 
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3.2 Fermions on the square lattice (U > At, |/ie| < f7 — 4^). 

Introducing the next order in the truncated hamiltonian, we obtain the antifer- 
romagnetic Ising model, lor which fi plays the role ol the magnetic field (see (2.25)). 
The phase diagram ol the truncated hamiltonian H^^^ consists ol three distinct regions, 
corresponding to the empty state, the two chessboard states, and the lull state. As 
before it is possible to bound the contribution ol the higher orders and thus to obtain 
part ol the phase diagram ol the effective hamiltonian, see fig. 1. The two separa- 
tion intervals (in thick lines) between the domains have length proportional to and 
make place lor domains corresponding to ground states with intermediate densities 
and higher periods. 

Let us remark that at this order the phase diagram is independent ol the magnetic 
fiuxes. 



— ] ^^ 



order 



or 



der 1 



1 ^ ' i k ^ ► 

U U 

Figure 1: Fermions on a square lattice. Phase diagram obtained at order U° and U~^. The dark 
domains represent values of /j, where no conclusion is obtained at this order. 

Since on the square lattice the expansion ol Hefj has only odd powers ol 
[/~^, the next order is U~^. Furthermore using the symmetry property ol the phase 
diagram we restrict the discussion to the case < 0. 

a) Study ofH^^K 

We first write the truncated hamiltonian H^^^ as a sum over potentials defined 
on 3 X 3 blocks (see (A. 3)). In Appendix A we give the zero potential Kb such that 
H'^ = Hb + Kb is an m-potential Irom which one obtains the following results valid 
il U is large enough. 



Let 



1 

/^i = -2- + (--cos(^) — 

_ f -2|- + (3 - cos(^)^ ilcos(^>0 
~ \ -2^ + (3 + f cos(^)^ ilcos(^<0 



1^3 



-2|^ + (3 + 3cos(^)^ ilcos(^>0 
fi2 if cos ^ < 



^2 fi 

/^4 = -2- + (y + 3cos(^)— (3.11) 
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then 

1. For /i < /ii, s_ is the unique ground state. 

2. For /ii < /i < /i2, H'^ is minimum for the configurations of fig. 2, and those 

obtained by rotations or refiections. 



o o o o • o 

o • o o o o 

o o o • o o 

Fi gure 2: Configurations defining Si 
3. For /i2 < /W < /W3, H'^ is minimum for the configurations of fig. 3. 



000 000 •o« 
o«o • o • 000 

000 000 o«o 

Figure 3: Configurations defining ^2 
4. For /is < /i < /i4, H'^ is minimum for the configurations of fig. 4. 



o o • o • o 

o • o •00 
• 00 o o • 

Figure 4: Configurations defining ^3 

5. For /U4 < /« < 0, H'^ is minimum for the configurations shown in fig. 5, and thus 
{■slfj^sli,} are the unique ground states. 



From the above discussion we conclude that the only configurations which yield 
the minimum energy for each block B are those represented on fig. 6, together with 
the empty and chessboard configurations, and those obtained by ion-hole inversion, 
and Proposition 3.1 is established. 
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o • o 
• o • 

o • o 



• o • 

o • o 

• o • 



Figure 5: Configurations defining the chessboard configurations 
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o 


o 


• 



Sl,p=\ S2,p=\ S3,p=l 



Figure 6: Ground states appearing in the phase diagram for fermions on square lattice; the occupied 
sites are black, while the empty ones are white. 

Proposition 3.1 ForU/t large enough, the ground states of the truncated hamiltonian 
at order 3 are the following Sj with density pj : 



QO < p < pi 




<S_ = {s. 


-} 




P- = 




1 

5 


Pl < p < P2 


Si : 


= {s?}, a = 


1,.. 


.,10 


pl = 


P2 < P < P3 




= {■Sal, « = 


= 1,- 


..,8 


P2 = 


1 

4 


P3 < p < P4 


Ss 


= {53}, « = 


= 1,- 


..,6 


P3 = 


1 
3 


P4 < p < 






slb} 




Pcb = 


1 




Note that S2 does not appear for cj) 

The phase diagram being symmetric with respect to // = 0, we have thus obtained 
the phase diagram for the truncated hamiltonian H^^^ (fig. 7). 

h) Ground states for Hpf f . 

We shall now show that the phase diagram for the effective hamiltonian is the 
same as the one for the truncated hamiltonian, except for domains of width of the 
order t^U~^ along the separation lines. Let us first make the following observation: 
for a chemical potential p in /ij+i [ the energy of any configuration s such that 
s\b ^ s^\b will satisfy the inequality 

H'^{s) - H'^{s^) > T > {) (3.12) 

where r is a function of p^ U and $. 
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Figure 7: Fermions on a square lattice. Pliase diagram for the hamiltonian truncated at order 3. 

In the following we consider the case fi G]/W3,/Wc6[ for which the six ground state 
configurations are given by S3 and to simplify the notation we write s" instead of S3. 
Let s be a configuration and Y a connected set of sites. We say that Y is ^^correcf 
if there exists some s" G S3 such that s\y = s"\y] moreover Y is said ^^a-correcf if 
there is a unique s" satisfying this condition. 

We then introduce the following notations : 

is the 3x3 block centered in x, 

Cp{x) is the p X p block centered on x if p is odd, centered on x + (|, |) if p is even. 
The integer p is the lowest common multiple of pi,P2, the periods of s in the 
directions 1 and 2. For example, p = 3 for S3. Note that if Cp{x) is correct it is 
necessarily a-correct with a = a[x). 

and for any configuration s 
r(s) = {x G A; Bx not-correct} 

B{s) = {Bx C A; Bx not-correct} \B{s) \ = \T{s)\ 



It follows that Cp{x) is not-correct if and only if Cp{x) C r(s). Therefore for any 
y ^ r(s), the condition Cp{x) 3 y implies that Cp{x) is a-correct and all Cp containing 
y will be correct with the same a[x). 

Let us decompose the complementary set r^(s) = A\r(s) into maximal connected 
subsets 




Cp(x)DB 




(3.13) 
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Then each Ik is a-correct, with a = a[k). 

With the definition of and (3.4), we have : 

BeB{s) 

+ E - (3.14) 

xeA 

Using then the inequahty (3.12) for the first sum and the expansion of R^^^ in terms 
of oriented paths, we obtain : 



1 /''+^ 

\Heffis) - i?e//(3«)| > r\Bis)\ " E E TT^ 

" xeAq>5 ^ 



(3.15) 



where because of (3.7) we have replaced s" by s"^^) (we set a[x) = a if Cp[x) is not 
correct). 

All circuits io such that io D Cp{x) ^ (f>^ io D T[s) = ^, are a(a;)-correct, and thus 



\Heffis) - H^n{s")\ > rms)\ - 4 E E ^ 

P xeAq>5 '■^ 



CO c ^ pupp L0\ 



(3.16) 



tjjr[Cp{x)i^4>, c^nr(s)7^(/) 
The function guj{s) is bounded by 2N[uj) • Then 

1 /''+^ 

\Heffis)-H,jfisn\>r\Bis)\--Y.-— 

P g>5 ^ 



|supp n Cp(x)| 



(3.17) 



c.efi,+i(.) iV(cu)|suppcu| 

The last sum is bounded by p^|supp c<;|, thus 

/g+i 

\Heffis)-H,jfisn\>r\Bis)\-j:-^ J: 1 (3.18) 

g>5 uj={xi,...,Xq^i) 
xieT{s) 

Let = I ^c'p{x)jB Cp{x)\^ [K = 9 in the case considered; K = 25 for 1S2, K = 49 for 
Si); then 

l^e//(3) - i?e//(3«)| > |B(3)| (^T - 9 (1)^(1 - f)-'^ (3.19) 

Therefore, for fi such that 

r>9(|)%(l-|)- (3.20) 

the set S3 is the set of ground states for the effective hamiltonian. In Appendix A we 
have shown the following result 
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Proposition 3.2 IfU/t > 387 then there existes^Ss such that for fi gJ/Ws + Es, /Wc6 — (^si 
the ground states of the effective hamiltonian are S3. 

For other intervals we proceed in a similar manner. 

3.3 Bosons on the square lattice 

From the expansion of the effective hamiltonian (2.28), one observes that at order 
1 the results for hard-core bosons are the same as for fermions. At order 3 the 

truncated hamiltonian for hard-core bosons and fermions are identical if (pp = 7r/2, 

(3) (3) 
moreover Hi (s; 6) = Hr ■ (s;7r — (h) if sp = 1, for all P. 

bosons^ ' ^ ■' lermions^ ' ^ ■' ^ ' 

We can write (2.28) in the form 
Hi^\s) = cte + J2Hp{s) + ^ E -^^-^y (3-21) 

P \x-y\ = 2 

Following the same method as for the fermions, one first shows that for U large 
enough and /« < 0, the ground states of H^^^ cannot contain bonds < x^y > with 

-1-2 -1-4: 

Sx = Sy = 1. Then one observes that, with jj, = —-jj -\- 6 



2U ' C/3 

00 



^p(oo) = -|7 + ^[4^-3] (3.22) 

i7p ( 82 ) = + ^[2^-2 cos (3.23) 

^p(oS) = -|7 + ^[ll-4cos<^] (3.24) 

Therefore for any configuration which does not contain bonds < x^y > with s^ = Sy = 
1, we have 

i7(3)(5; S, ^) = H^^\s; S - cos^,^ = + cte((^) 



<Uons(^;'^-c°«<^'<^=i) + cte (3.25) 



(where cte(^) is independent of s and 6) and thus the phase diagram for bosons with 
hard-core is directly obtained from the phase diagram of fermions with (f> = (fig. 
8). 



3.4 Fermions on the triangular lattice 

As mentioned in the introduction the triangular lattice is especially interesting 
since it is not bipartite and the system does not have the usual particle-hole symmetry. 

From the expansion (2.26) for the effective hamiltonian, one sees immediately that 
at the order the truncated hamiltonian H^^^ is the frustrated antiferromagnetic 
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Figure 8: Hard-core bosons on a square lattice. Pliase diagram of the hamiltonian truncated at 
order 3. 



Ising model with magnetic field Going to next order, which is U~'^ in this case, 
we will see that the frustration disappears. 
At the order U~'^ we have : 



H,f f is) = H^^Hs) + R^^Hs) + cte 



(3.26) 



Ha{s) = - (y^ + ^ cos <P^^ E^^ + E ^^^y + COS <P^s^. 

\ / xE.A <a;,j/>CA 

Computing Ha[s) for the eight configurations on the triangle A, we have the following 
result; let us define 



1^2 



-6— cos (pA 



(3.27) 
(3.28) 
(3.29) 



/«3 = 3-- + 3— COS^A 
U U'^ 

then the ground states of are given by (see fig. 9) 

H < Hi 71 = /?_ = 

Hi < ji < Ii2 Ts = {s^}, a = 1, 2, 3, P^ = \ 

2 

H2 < H < H3 Ts = = 1,2,3, /05 = - 

fi> fi3 T+ = p+ = 1 

where Sg is a configuration with one of the three sublattices occupied, and Sg is a 
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configuration witli two sublattices occupied. 
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o«o«o«o«o o«oo«oo«o o««ooo««o 

oooooooo oo«oo«oo ooo««ooo 

o«o«o«o«o o«oo«oo«o o««ooo««o 



p — \ T5, p — ^ %, P 



o«oo«oo«o o««oo««oo 000000000 

o««o««o« oo««oo«« •••••••• 

o«oo«oo«o ©©oo^^oo* 000000000 

oo«oo«oo o««oo««o •••••••• 

• •o««o««o ©oo^^oo** 000000000 

oo«oo«oo ©©oo^^oo •••••••• 

o«oo«oo«o oo««oo««o 000000000 

o««o««o« voovvoo* •••••••• 

o«oo«oo«o o««oo««oo 000000000 



4 _ 1 _ 1 

9 ■'St P — 2 ■'9t P — 2 



Figure 9: Ground states appearing in the phase diagram for fermions on triangular lattice; the 
occupied sites are black, while the empty ones are white. 



Tlie pliase diagram for tlie truncated liamiitonian H^^^ is represented on figure 10, 
wliicli sliows explicitely tliat tliere is no particle-liole symmetry. On tlie boundary tlie 
ground states are infinitely degenerate. Following the method described in section 3.2 
one can then prove that for U /t sufficiently large the phase diagram for the effective 
hamiltonian remains the same, except for domains with width of order /U'^ around 
the boundaries. 

Going to order 3 the problem of the ground states of the truncated hamiltonian 
becomes more difficult. In the Appendix we give zero-potentials from which one 
can deduce the ground states of H^^^ for some domains of the (//, ^)-plane, but we 
could not find a rigorous proof for all values [p^(f>). Let us note however that except 
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for the configurations '7^, '7^, all the configurations 7^, a = f, . . . ,9 (fig. 9) do appear 
in the domains where the ground states have been rigourously established. 
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Figure fO: Fermions on a triangular lattice. Phase diagram of the hamiltonian truncated at order 
2. 

The results discussed in the Appendix are summarized as follows. Let us define 
the following functions of (f> 



fi7 



-3— + 3 cos 




-3— + 3 cos 



-6 cos <p— + 

0"^ 



-6 cos <p— + 



-6 cos (f)— + 



[/2 



(-f + fcos2^ 
(-f + 9cos2(? 

9 9 

cos 2(i> 

i6 4 ^ 



(-1 + 1 cos 2^ 
(_| + |cos2^ 
(-§ + 9 cos 2^ 



) U3 

' (763 



if 
if 



COS 

5_ 
12 

if cos 



> 



< cos ; 



< 



^ 12 



1 



4 

if cos ! 



cos . 



< 



5_ 
12 



> 



'(73 



if I < COS ; 
if < cos 26 < 



Yn if COS 2(/) < 



< 



(3.30) 
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Figure 11: Fermions on a triangular lattice. Phase diagram for the hamiltonian truncated at order 
3, for chemical potentials corresponding to low densities. 



With the configurations 7^, a = 1, . . . , 9 represented on fig. 9, the ground states 
of H^^^ are given around // = — by 





P- = 


-- 












Pi = 


1 

7 


Hi < n 


< 


P2 






P2 = 


1 

5 


< P 


< 


{ 


P3 














Pa 




P3 = 


1 

4 


< P 


< 


Pa 


and 




P4 = 


1 

4 


P3 < P 


< 


Pa 


and 




P5 = 


1 

3 


Hi < n 


< 


P5 





1 

2 



while for values of the chemical potential around 0, 

/i6 if 4 < cos 2^ < I 



p7 

Tg, P9 



P5 < P < ] -r / 1 

[ p-r il cos Z(p < ^ 

j ps < p < pr if cos 2(f) > 
\ pe < p < P7 if 4 < cos 2 
P7 < p < P7 and cos 2(f> > 
P7 < p < P7 and cos 2(f> < 



^ 12 



where p^ is defined from u-, as follows : the coefficient of tpj is identical, while the ones 



Pj 



C/2 



of jj and jj^ are opposite. 



C/3 
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The phase diagrams are represented on figures 10, 11 and 12. 




Figure 12: Fermions on a triangular lattice. Phase diagram for the hamiltonian truncated at order 
3, for chemical potentials corresponding to densities close to one half. 



3.5 Bosons on triangular lattice 

From the expansion (2.29) one observes that at the order U~^^ the truncated 
hamiltonian for hard-core bosons is the same as for fermions (antiferromagnetic Ising 
model with field). Going to the next order, i.e. we have : 

from which one concludes immediately that the ground states of H^^^ are the following: 



/i < -3— - 3cos^— 

-3-- - 3cos < /i < T5 
U U'^ 

- 

< /i < 3-- + 3cos ^— T5 
/i > 3— + 3cos^— 



T_={3_}, p. 

{■55} « = 1,2,3, p^ : 

{s^} a = 1,2,3, pr, -- 

^+ = {■5+}, P+ 
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Figure 13: Hard-core bosons on a triangular lattice. Phase diagram for the hamiltonian truncated 
at order 2. 

Then except for domains of order ^ the phase diagram of the effective hamilto- 
nian is the one presented on fig. 13. 

In the Appendix we then show that it is possible to find zero-potentials such that 
the truncated hamiltonian at order 3 is an m-potential. We then obtain the following 
result. Let 

|-3cos2<^)^ (3.32) 
;^-3cos2<^)^ (3.33) 
3cos2(^)^ (3.34) 



jll = 


4- 


- 3 cos (f. 


' h 


1^2 = 


-4- 


- 3 cos (f. 


' h 


1^3 = 


-4- 


- 3 cos (f. 


' h 

U2 ^ 


fi4 = 


9 

'-5- 


3 cos 2(1 




1^5 = 


|(-2 


— 3 cos 

— 1 cos 


2<^)^ 



4 ^'[/3 



if COS 2(^ > 
if cos 26 <0 



(3.35) 
(3.36) 



Then the description of the phase diagram of the truncated hamiltonian at order 
3 is as follows 





p- 


= 


fi < fil 




pl 


_ 1 
~ 7 


jll < jl < 1^2 


^1/45 


p 


_ 1 
~ 4 


jJ,2 < P < P3 




P5 


_ 1 

~ 3 


jl3 < jl < jl4 




P6 


_ 2 
~ 5 


jl4 < jl < jl5 
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Figure 14: Hard-core bosons on a triangular lattice. Phase diagram of the hamiltonian truncated 
at order 3. 

7^, /08 = | /«5</«<0 and cos 2^ < 

7^, /J9 = | /W5</W<0 and cos 2^ > 

where T^j^ is the set of all configurations with density | and the constraint that 
with respect to a lattice direction : each two lines is empty, while the others have 
alternatively one ion and one empty site (log I'T^'^^I = (9(|A|2)). This set is probably 
not a set of ground states for the effective hamiltonian. 

The phase diagram is symmetric around // = 0; therefore the separation lines for 
/i > are derived immediately from /ii, . . . , /i5. 

4 Flux phase problem 

4.1 Minimum flux for a given conflguration of ions 

In this section we consider the following problem : given a configuration s of ions, 
what are the fiuxes across the elementary cells C which minimize the effective 

hamiltonian ? As in the previous sections, A is either the square or the triangular 
lattice, and the cells C are either plaquettes P or triangles A. The proofs are based 
on the expansion, thus the chemical potential for the electrons must lie in the gap 
and U has to be sufficiently large. We take again /ie = and write //j- = 

The case of bosons was treated in section 2.4. From now on we consider only 
fermions systems. 

Theorem 4.1 1. For the configurations s = or the effective hamiltonian is 
independent of the fluxes. 
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2. Let s ^ satisfying the condition that for any site x there exists a site 

y with \x — y\ < d[s) < oo and s^ = —Sy. Then there exists a positive number 

Uq (depending on s) such that for U > Uq, the minimizing fluxes are given by 
(^mm^^^ ^ ^^^^^ ^f^^^^ 

. 7 / \ I TT (f ^ sites of P are occupied / , \ 

square lattice : (pp{^) = i 7 ■ l'*-!) 

[ otherwise 

. - \ IT if or 1 site of A is occupied , , 

triangular lattice : (pA = i . r \ ■ ? l"*-^) 

10 if 2 or 3 sites of A are occupied 



The detailed proof of theorem 4.1 can be found in Appendix B. Here let us outline 
the main ideas. 

1. If s = s+ or s_, 0(s) is empty and the effective hamiltonian does not depend 
on the fluxes. This is not surprising, because we took the chemical potential of the 
electrons to be in the gap, and therefore the number of electrons is either or |A|; 
in the first case there is no electron to feel the magnetic field, and in the second case 
the electrons cannot hop because there is no free lattice sites; then the phases of the 
hopping coefficients have no infiuence. 

2. Given s ^ s_} defined on the square lattice, let V denote the set of plaque- 
ttes which are neither completely occupied nor empty. From (2.17) the contribution 
involving fiuxes in the truncated hamiltonian H^^\s) is 

2 ^ . _ , 2 2 t* ^ 2 , , 

^P- TTTT:^ TT = -77T >. cos^p- ^ (4.3) 



— Ti 

8 [/3 



■1 



COS . 



m 



l)!(3-m)! 8[/3^''°"^''(m-l)!(3 



m 



where 



Lpp 



Dvr; 



-1^ 



Sp 



(4.4) 



Therefore the fiuxes which minimize i7(3)f5) are <f)p = |(sp + l)7r if P G P and (f)p 
arbitrary if P ^ V. 

We then consider the truncated hamiltonian H^^\s). Let TZ denote the set 
of rectangles (two plaquettes) which are neither completely filled nor empty. We 
decompose 7^ = 7^l U 7^2 with 7^l = {Pi = Po U P'; Pq ^ V,P' G V} and 
7^2 = {P2 = P' U P"; P' G V,P" G V}. The contribution involving fiuxes in 
H(^Hs) is 



2 t^ 
SIP 



E 



1 



COS (ppi 



2 f_ ' 

(m'- l)!(3-m')! ^ U^^\ 



C2 COS tpp^ + — — COs((^p/ + tpp^ 



411 



1 

+ — cos (^p/- 



I L _c' 

(m'- l)!(3-m')! IP ^ 



+ 



(m'- l)!(3-m')! [/2 ^ 



+ 



-S^y COS((^p/ + (^p//) 



(4.5) 
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where xy = P' D P" and we used the fact that spispn cos(^p/ + cf>pii) = cos((^p/ + (,cp"). 
We note also that I C\ I , I C-^ I , I C-^ I QjrG smaller than 3 and C — 2 > 0. 

Therefore for U > At the first term is minimum for (ppi = |(sp' + l)7r as before 
and (fiPg = 0. On the other hand the second term is of the form 

a cos (fi' -\- l3 cos (fi" + e-rrr cos((^' + (^") 

where a, /3 > |(1 — 3^) and |e| = 1. Therefore for U > At this term is minimum 
for (fipi = (fipn = TT, as was obtained already at the previous order. Finally the fiuxes 
on those plaquettes P such that P fl 7^ = remain arbitrary. We then proceed by 
induction to determine the value of (pp for the remaining plaquettes. 

4.2 Flux phase problem 

We first discuss the fiux phase problem for the truncated hamiltonian at order 
3. For a given configuration s, the fiuxes {^(^(s)} which minimize H^^\s^ {^c}) are 
given by (4.1) and (4.2) when C is not filled or empty, and arbitrary otherwise. 

Introducing appropriate zero-potentials (see Appendix C) we may write 
{(f>c{-s)}) as an m-potential 

H('\s; fi; {Ms)}) = E H'b{s; p; {Ms)}) (4-6) 

B 

and one can then find its phase diagram, and thus the fiux phase problem is solved 
for the truncated hamiltonian at order 3. 

Theorem 4.2 For U sufficiently large, the configurations of ions and the set of fluxes 
which minimize the truncated hamiltonian at order 3 are the following (fig. 15 and 
16). 

i) for the square lattice 



jl < jll = 


t^ 1 t* 

~^[/ ~ 2[/3 




P- = 


0, , 




is arbitrary 


jll < jl < jl3 = 


t^ t^ 

-2— + 2 — 

U 


Si, 


Pl = 


1 

5' ' 


j)p 


is or arbitrary 


jl3 < jl < jlcb = - 


t^ 13 t* 

-2— + 

U 2 


Ss, 


P3 = 


1 

3' ' 


j)p 


is or TT 




jlcb < P <0 




Pcb = 


1 

2' 


(j)p 


is TT 



Using the particle-hole .symmetry the phase diagram for fi > is obtained by 
replacing s„ by —Sa, without changing the fluxes; 
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a) for the triangular lattice, 

t^ t^ 9 

fi < fii = -3— - 3— - -— T_, p_ = 0, <pA is arbitrary 
U (7^ 4 ty-^ 

t^ t^ bt^ ^ 1 7 

jii < ji < Ii2 = -377 - 3— + -— Ti, /?! = -, «s TT or arbitrary 

U U 4: U'^ I 

t^ t^ Q9t^ ^ 1 ^ 

< IJ. < IJ.3 = -3— - 3— + — — T3, /J3 = -, tS TT 

15 ^ 1 7 . 

iJ-^ 5 
jiQ < ji <{) /jg = -, ^A «s or TT 

(4.7) 

T/ie phase diagram for fi > is obtained by replacing ^„ 6j/ — — 

^„ - TT. 

Remark : for the square lattice, the configurations S2 do not appear except at 
the value // = /is which is infinitely degenerate. In the case of the triangular lattice, 
the configurations 7^, 7^, and Tg do not appear. 



S- Si 1S3 Scb Si iS+ 
\ \ \ \ \ ^ ^ 

1^1 1^3 f'ch —f'ch —1^3 —f'l 

Figure 15: Flux phase problem on square lattice. Phase diagram of the hamiltonian truncated at 
order 3. 



T_ Ti % T5 Te % 
\ \ \ \ \ ^ ^ 

jll jl2 /«3 1^5 1^6 

Figure 16: Phase diagram of the flux phase problem on triangular lattice. 

Following the ideas developped for the case of constant fiuxes, it is possible to 
define subintervals where a suitable Peierls condition can be obtained for 

H^^\ and to prove that the phase diagram for the effective hamiltonian is the same as 
the phase diagram for H^^^ up to domains of width of the order ^Ja (or in the case 
of the square lattice) around the boundary points. 
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5 Conclusion 



The phase diagram at finite temperature has been obtained in [12, 16], only for 
the square lattice and around the symmetry point. There the two chessboard ground 
states are related by a particle-hole or translational symmetry and therefore a Peierls 
argument can be applied. For other ground states this is not the case and Pirogov- 
Sinai theory has to be used^. Recent developments of this theory allowing general 
long-range (and also quantum) interactions may be found in [4] and [3]. Using these 
results it is expected that the phase diagram is stable at low temperatures. 

The present work is limited to the large U limit. For small U the phase diagram 
changes drastically. This is the case in one and two dimensions where a very rich 
structure has been uncovered [9, 5, 27]. We expect that values different from or tt 
will appear for the optimal fiux as U is lowered. 

In the case of hard-core bosons when [/ = 0, we get the quantum XY model. 
This model is known to display off-diagonal long-range order (ODLRO) at half-filling 
in the ground state [13]. On the other hand as soon as U > the ground state 
configuration for the nuclei is a chessboard (on a square lattice) [22]. For large enough 
U it seems clear that the bosons will form a charge density wave of period \/2- An 
interesting question is whether the critical value Uc for the disappearance of ODLRO 
isUc = or Uc > 0. 

Finally we wish to point out the link between the Falicov-Kimball model and 
Hartree-Fock theory for the Hubbard model. If one minimizes the ground state (or 
free energy) for the repulsive Hubbard model, only over the set of Hartree-Fock states, 
then one finds that the mathematical structure of this variational problem is very 
similar to the one for the Hubbard model. This fact has been exploited in great detail 
recently [2]. The results obtained here on the ground states for rational filling factors 
of the form ^ may shed some light on the Hartree-Fock states that one should look 
for in the Hubbard model. 

We hope to come back to these questions in the future. 
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A m— potentials 



A.l Square lattice 

Following the method of [8] using plaquettes and bonds < xyz >, it is straight- 
forward to show that for U > 4t, one has the following property : for // < /ii, s_ is the 
only ground state of H^^\ and for > // > /ic6, <Sc6 = {-slbT^lb} only ground 

states of H^^\ where 

1 

/^i = -2- + —{- - cos <P) 
P'cb = -2- + Y + 3 cos (f)) 

To obtain the ground states of H^^^ in the remaining interval // G [/ii,/ic6], we 
first note that for U > At and fi < ficb the ground state configurations cannot contain 
bonds < x^y > with = Sy = 1. 

Introducing the parameter 6 defined by 

we write the truncated hamiltonian in the form 

H^'Hs) = ^^ ^ l + J^J2HBis;^) + cte (A.2) 

^ <x,y> BCA 

where B are 3x3 blocks, centered on x(B), and 



Hb{s; fi) = -S{s^^B) + i) - + 2 cos (f)) ^ s^Sy 

° <x,y>CB 

+:^(4-cos(^) ■^xSy + ^ + :^cos(^ ^ (1 + 5sp) (A.3) 

\x-y\=V2 \x-y\ = 2 

As explained in the introduction we then introduce zero potentials Kb such that 
H'g = Hb + Kb is minimum if and only if the configurations on B are either S\ , S2 or 
1S3 depending on the value of // (i.e. (5). Numbering the nine sites in a block B from 
left to right, starting down and ending up, let 

KB = i2^^k'i^ (A-4) 

8 = 1 
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where the are given in Table 1 and 

kg^ = -52 - S4 - Se - Sg + 4S5 (A. 5) 

= -si - S3 - sr - sg + 4s5 (A. 6) 

k^B = -^1^3 + 2S4S6 - S7S9 - S1S7 + 2s2S8 - S3S9 (A. 7) 

kg^ = —S2S4 — S2SQ — S4SS — SqSs + S1S5 + S3S5 + S5S7 + S5S9 (^-8) 

These potentials k^^^ are invariant by reflection and rotation. 

Doing explicit calculations, we find that the only configurations which minimize 
H'g are the following : 

If cos ^ > 0, then for 

I — cos (/) < S < 3 — cos (/), it is {s^} 
3 — cos (f> < 6 < 3 + 3 cos ^, it is } 
3 + 3 cos ^ < (5 < ^ + 3 cos ^, it is {s^} 

If cos ^ < 0, then for 

I — cos ^ < (5 < 3 + I cos <f>, it is {s^} 
3 + I cos ^ < (5 < Y + 3 cos ^, it is {s^} 

which establishes the Proposition 3.1. 

We then have to find intervals + e^, fij+i — Sj[ such that 

H'b{s) - H'sis]) > T >0 (A.9) 

with 

and K the constant defined in section 3.2. 

Explicit computation for the case fi G]/W3,/Wc6[ discussed in sec. 3.2, yields for 
cos ^ > 

.3 = 27.2-(l-|)"'^ (A.U) 

<,, = 27.2»(l-|)"^ (A.12) 

and the condition 63 -\- 64 < jicb — /"s = |(3 + cos will be satisfied ii U > 387. 

In the case of the bosons, as discussed in section 3.3, we can use the results for 

— 2L 
2 ■ 



'7 

the fermions with constant external magnetic fiux 



A. 2 Triangular lattice 
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«! a2 C(3 Ct4 



Si , COS ^ > 


1 . rrit; m . OS 
j_ 1 ' ' ' ~ 1 ^ (J 

40 ~'~ 120 ~'~ 15 


1 . rrit; m . S 

J_ 1 ' ' ' ~ 1 u 

480 ~'~ 40 ~'~ 15 


1 rrit; m S 

480 240 240 


I 

~32 


Si , cos ^ < 


11 3cos(/) 1 25 
240 40 ~'~ 15 


1 1 COS (/> 1 S 
120 ~'~ 15 ~'~ 15 





1 COS (f) 

24 24 


S2 , cos ^ > 


16 ' 8 


1 1 cos (f> 1 S 
24 ~'~ 48 ~'~ 16 


1 COS (f> 1 5 
96 32 ~'~ 96 


1 

16 


1S3, cos ^ > 


5 
8 


1 1 cos (f> 1 S 
96 ~'~ 48 ~'~ 16 


1 cos (f> 1 5 
24 32 ~'~ 96 


1 

32 


1S3, cos ^ < 


5 
8 


5 cos (f> 1 5 
96 96 ~'~ 12 


11 3cos(/) 1 5 
480 80 ~'~ 240 


9 7cos(/) 1 5 
160 160 ~'~ 120 



Table 1: Coefficients for the definition of the zero potential Kb- This is for the truncated hamiltonian 
up to order 3 in the case of fermions on a square lattice, fi = — — + jjsS. 



Fermions 

The 3x3 blocks of the square lattice are replaced by hexagons H : the hexagon 
centered in x is defined as{j/GA:|a; — j/|<l}. The truncated hamiltonian may be 
written as 

A y a;GA <a;,j/>CA 

t"^ COS (h 3 cos (h t"^ \ ^-^ , , . . X 

-^71^1:^. + ^^^-) E «»(-«) + ct.e (A.13) 

where 

7 + 5cos2^ , , 4 — cos 2^ , , 

^H{s) = — 2^ S^Sy + — 2^ S^Sy 

<x,y>CH (x,y)CH 

\x-y\=\/?> 

1 ^-^ 5 COS S 

+ 7 E ^xSy + — — ^ y sp (A.14) 

(x,y)cH PCH 
\x-y\ = 2 

As before we introduce zero potentials k^^] numbering the sites starting from the 
center, i.e. 1 is the center and 2-7 are the external sites, we define 

k^H^ = (^si-ys, (A.15) 

8 = 2 

7 6 

^H^ =^SiS, -yS,S,j^l - S2S7 (A. 16) 

i=2 8 = 2 

5 5 

^H^ = Si{yS,S, + 2 + S6S2 + S7S3) - yS,S, + iS, + 2 - SeS7S2 - S7S2S3 (A. 17) 

8=2 8=2 

We set as before H'^j = Hh + Kh with Kh a linear combination of k^^ and 

f2) 

k\j such that H'fj{s; ji) is an m-potential. Hereafter we define domains in the phase 
diagram, for which we can find an m-potential Kh = J2i (^ikjj with ai given in Table 
2. Unfortunately these domains do not cover the whole phase diagram, and thus there 
are still regions for which our results are not proven. 
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I)_ : domain of 71 



T>i : domain of Ti 
T>3 : domain of T3 

: {{f^,<P);f^cb < fi < -f }, with 




-3|^ + 3cos(^^ + (f + 9cos2(^)^ ifcos2(^>0 
-3|- + 3cos(^^ + (f + 1 cos2(^)^ ifcos2(^<0 



(1)5 is a part of the domain of 7^) 

Df^ : {(/i,(^);cos2(^ < ^ and - < // < -6 cos + (-f + 9 cos 2(^)^} 

I^g"^^ : {{fi, <f>); cos 2<f> > < fi < fis, and also : 

^ < -6 cos ^f, + (-^ + '-^)^ and < -6 cos + {% - '-^)^} 
(1^5^^ and I^g"^^ are also parts of the domain of 7^) 

De : domain of Tq, under the restriction // < —6 cos (f)-^ + (| ~ ^^^y^)^ 

Di : domain of Ti, with // < —6 cos (f)-^ + (| ~ ^^^^) 

I)g"^^ : domain of Tg, with fi > —6 cos + ( — | + 3 cos 2<f>)jjj 

V^s^ : domain of Tg, with // < -6 cos (f)^ - ^^^^ 

T>g : domain of Tg 



«! a2 C(3 



D_ and Di 


45 19cos2(/) 5' 
224 224 7 










1 cos2(/) 5' 
8 32 8 










1 5cos2(/) 5' 
96 32 9 










3 3cos2(/) 5 
128 32 8 


11 9cos2(/) 1 5 
64 32 ~'~ 12 


17 cos2(/) 1 5 
384 32 ~'~ 72 




1 3cos2(/) 55 
224 112 42 


47 33cos2(/) 1 25 
224 224 ~'~ 21 


17 1 9cos2(/) 1 5 
224 ~'~ 112 ~'~ 42 




5 
8 


33 87cos2(/) 1 5 
160 160 ~'~ 20 


3 9cos2(/) 1 5 
80 40 ~'~ 40 




1 cos2(?!) 6 
32 16 9 


11 1 cos2(;i 1 5 
32 ~'~ 32 ~'~ 18 


1 cos2(?!) 1 5 
32 16 ~'~ 18 


^8 


5 
8 


9 3 cos 2?^ 1 5 
32 32 ~'~ 12 


5 

24 


V'i^ and Dg 


5 
8 


9 3cos2(/) 
32 32 






Table 2: Coefficients for the m-potential. We defined 6 s.t. fi = —6 cos (f)jj2 + 6-^ and 6' s.t. 

— _3*i _L 3 cos i/it^ _1_ t* CI 

P — jj -I- JJ2 "T C/3 • 
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«! a2 C(3 





47 5cos2(/) (5 
112 28 7 








/i2 < /« < /«2 


11 cos2(/) (5 
32 8 8 








/ / 


11 cos2(/) 5 
48 12 9 


U 


U 




A* 
8 


11 cos2(/) 1 /i 
16 8 ~'~ 12 


cos2(/) 1 /i 
24 ~'~ 72 


fi4 < fi < (—2 — 3 cos 2<f>)jjj 


A* 
8 


43 9cos2(/) 1 /i 
80 40 ~'~ 20 


1 1 3cos2(/) 1 /i 
20 ~'~ 40 ~'~ 40 


(-2 - 3cos2(^)^ < /« < 


A* 
8 


7 3cos2(/) 
16 8 






Table 3: Coefficients for the m-potential of bosons on triangular lattice. We defined 6 such that 
fi = ^ '^^2*^* — ^ Ti^^- T^^^ * = 2, 3, 4 were previously defined. 



Bosons 

The construction of the appropriate m-potential is given in Table 3. 

B Proof of the theorem 4.1 



We will show that for any configuration {(f>c + ^c} with 6c ^ for at least one 
cell C we have 

H,jj{s; fi; {^c + Sc}) - H,jj{s; fi; {^c}) > (B.l) 

Recall that Hefj is given by (2.16) where 0(s) is the set of circuits io which visit at 
least one empty and one occupied site. We introduce the set of circuits 0(s) such that 
id G 0(s) and there exists a cell C^o which is enclosed by io but not by any shorter path. 
We remark that for any io G 0(s) the cell C^o is unique and N[ijj) = 1 (Eq. (2.5)). 
We also denote by kc{io) the winding number of u around the cell C . For u G ^{s), 

kc{L0) = ±1. 

We split the left hand side of (B.l) in two parts /i + /2, where /i comes from the 
circuits u G 0(s) and I2 from the circuits u G 0(s)\0(s). For /i we have 



cos <^ ^kc{to){<j)c + 



c 



cos 1^ kc{io)ct>c 



:b.21 



This term is decomposed further as /i = Ji + J2 with 



C^=C 



COS < ^ kc'{oj)<pc' + kc{oj)6c 

I C" 



COS < ^ kc'{oj)<pc' 
I C" 



:b.3) 
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and 

cos < ^ kc'{to){(i)c' + 



c^=c 



cos < ^ kc"{io)(f)c" + kc{to)8c 



:b.4) 



For I2 we have the same formula than (B.2) with the sum over u G 0(s) replaced 
by the sum over u G 0(s)\0(s). In what follows we show that Ji is bounded below 
by a strictly positive number and that J2 are of smaller order. 

Lower bound for Ji 

Since Y,c" kc'{io)(f>c' = or tt and kc'{io) = ±1 the difference of the two cosines in Ji 
is equal to 

cos 1^ A;c'(cu)(^c'| (cos (5c' - f) (B.5) 
Therefore 

Ji = J2'^c{s){l- cos 6c) (B.6) 

c 

with 

^cis)= E 77S^k(^)l(-ir^^^-'cos(^^a'H^C'| (B.7) 

c^=c 

One can then show that for any io G 0(s) 

(_1)™H-1cos|x;A;c'H'^c'| = 1 (B.8) 

and therefore ^cl'S) > 0. We give the details of the proof of (B.8) only for the square 
lattice, the case of the triangular lattice being similar. 

Ceo is the unique cell surrounded by cu, for which is a circuit of minimal length. 
Then necessarily io is contained in a rectangle with C^o in one corner, and the occu- 
pation of all sites in the rectangle is the same except for the site x (or maybe the two 
sites X and y) in the corner opposite to C^o (fig. 17). 



Therefore in the case of figure 17 (a), (pc' = for all C surrounded by cu, and 
( — 1)™('^) = —1, which implies (B.8). On the other hand, in the case of figure 17 (b), 
(f)C" = for all C ^ C\ surrounded by uj and (f)c\ = |(1 + SxSy)Tr. Thus if uj contains 
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Fi gure 17) Minimal circuit for an empty square; and/or Sy = —s. 

only one of the sites (x^y) (B.8) follows as previously , while if contains both sites 
X and J/, 



cos A;c'(tu)^c'| = s^Sj^cos|^(l + s^Sy)Tr^ 

Since m is either 1 or j — 1, j — 2, we have in all cases 
1 



-1. 



where j{C) is the length of the shortest path io G 0(s) surrounding C. Since (1 
cos Sc) > '^^c ^6 from (B.6) 



B.9) 



1 1 fJ^^^ 



c 



2 2J(c?)-2 U0(^')- 



c 



:b.io) 



Now we proceed to find upper bounds for I2 and J2. 



Upper bound for I2 



Using |5'c^(s)| < I for all cu, |sina| < |a|, and cos a — cos /3 = —2 sin sin ^y^, we 



may write 

c^efi(s)\fi(s) 

<- Y — — 

c^efi(s)\fi(s) 



Y.\kcH\ 



c 



:b.ii) 



where the second inequality comes from Cauchy-Schwartz. We bound Xlc l^c(^)| by 
7|c<;p, with 7 a constant depending on the lattice (square lattice : = j^] triangular 
lattice : 7 = |)- Then 

,0 tH 



<^ c^efi(s)\fi(s) 



:B.12) 
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We have |A;c(c<;)| < j^. Moreover since the circuits belong to 0(s)\0(s) and surround 
C their minimum length is j{C) + 1; therefore the number of circuits of length n 
enclosing C is less than 7'n^2;"~^, with the maximum coordination number z and the 
constant 7' depending only on the lattice. Thus we obtain 
„/ 



^\^\ C n>]{C) + l 



r 



n-1 



n>l 



u 



:b.i3) 



with the hypothesis of theorem 4.1 there exists jmax such that j{C) < jmax and 
therefore 



t P^^^ 
where the constant a does not depend on 



:b.i4) 



Upper bound for J2 

Proceeding as in the previous paragraph we write 



|J2|< E -^^2 



2lM- 



1 



E kc{io)8c 



1 



The first sum may be bounded as previously and 



:b.i5) 



:b.i6) 



We remark for later use that in the last sum the lengths of circuits are equal to 
j(C"), and since they surround C we have j(C) < j{C'). 



Proof of (B.l) 

Collecting the three bounds (B.IO), (B.14) and (B.16) we get 

Heff{s; fi; {(pc + ^c}) - Heff{s; fi; {<f>c}) = Ji + J2 + h 

c,c' c^efi(s) 

Ct^C" c^=C' 



c 



where 



1 1 



t 

2a- 



:b.i7) 
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Let K{C) be the number of occurences of C in the second sum; then K{C) < 
7"imax5 with 7" depending only on the lattice, and 

HC) fA^) f_b{C^J^., 
- ^ K{C) + 1 U^i^)-^ ^ ^ ^, \ K{C) + 1 U^i^)-^ ^ 

j(c?)<j(c?') 



6(C) 



1 



K{C") + 



+ E 

j(c?)<j(c?') 



HO , 

^HEll'^^l 



\ + 1 u' 



HC) 



\ K{C") + 1 u 



:B.181 



In writing this last inequality we used 



c^=c" 



t 2 



:b.i9) 



'{K{C) + 1){K{C') + I) 
which is true for U large enough, since j{C) < j{C'). Theorem 4.1 follows from (B.18). 



C Flux phase problem 



C.l Square lattice 

From (2.15) we obtain an expression similar to (A. 2) with 
Hb{s; fi) = -S{s^lB) + + E ^^^y 

\x-y\ = 2 



'^^ PCB 



1 - E ■'^^■'^y + ^■'^P I cos - - E ■'^^■'^y + 4 E ■'^^■'^y + 2 

(x,y)cP ) ^<x,y>CP {x,y)CP_ 

\x-y\=-^2 



(C.l) 



As before we consider the new potential H'^ = Hb + ^1^=1 (^i^^B with the coefficients 
ai defined in Table 4 and the zero-potentials k^g given by Eq. (A.5)-(A.8). One can 
check that H'g is an m-potential leading to the phase diagram of fig. 15. 
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«! a2 C(3 Ct4 



Si 


s 

8 


1 ^ 5 
16 ' 16 


7 ^ 5 
32 ' 96 


1 

32 


Ss 


3 _|_ 25 
80 ~'~ 15 


29 _|_ 5 
480 ~'~ 15 


1 

4 


1 

24 



Table 4: Zero-potentials Kb for the fermions on a square lattice, with minimal flux. 



C.2 Triangular lattice 

When /i is close to — — 3^ the ground state configurations consist of triangles 
containing zero or one ion; otherwise the energy increases as jj. 

The minimum fiux for these configurations is tt everywhere; since it is constant 
we can use the results of section 3.4 and thus the domains of 71, 7^, T3 and part of 
T5 are proven. 

For the remaining part we write the truncated hamiltonian as 



E - 



<x,y>CA 

J2 Hh{s]ij) 



3sa — cos E 

X^A 



HCA 



(C.2) 



with 



Hh{s; fi) 



— T s s - V 

"32 ^ ^^rr'^""'^" ^4 ^ 



\x-y\=V3 



(x,y)CH 
\x-y\=2 



1 



+ — E^o^^-Pp-" E -Sx-Sy + bsp 

PCH \ (x,y)cP , 



(C.3) 



where = Eacp ^a- 

Adding to Hh[s; //) the zero-potentials defined by (A.15)-(A.17) with the factors 
given in Table 5, we obtain an m-potential proving the phase diagram of fig. 16. 



«! a2 C(3 



% 


A* 
8 


7 

32 


+ 


A* 
12 


1 _|_ A* 
24 ~'~ 72 


% 


A* 
8 


15 

32 


+ 


A* 
20 


1 _|_ A* 
8 ~'~ 40 


% 


A* 
8 




23 
32 








Table 5: Zero-potentials Kb for the fermions on a triangular lattice, with minimal flux. 
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